Antiferromagnetic hedgehogs with superconducting cores 
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Excitations of the antiferromagnetic state that resemble antiferromagnetic hedgehogs at large dis- 
tances but are predominantly superconducting inside a core region are discussed within the context 
of Zhang's S0(5)-symmetry-based approach to the physics of high-temperature superconducting 
materials. Nonsingular, in contrast with their hedgehog cousins in pure antiferromagnetism, these 
texture excitations are what hedgehogs become when the antiferromagnetic order parameter is per- 
mitted to "escape" into superconducting directions. The structure of such excitations is determined 
in a simple setting, and a number of their experimental implications are examined. 
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I. INTRODUCTION 



In Zhang's S0(5)-symmetry-based approach to the 
physics of high-temperature superconducting materi- 
als the local state of the system at the spatial point 
r is characterized by the orientation of a five-component 
unit vector n(r). Orientations for which $^ =i ( n<1 ) 2 = 
are purely superconducting, the orientation of n in the 
(two-dimensional) 4-5 hyperplane determining the phase 
of the complex superconducting order parameter. Ori- 
entations for which Y^a=4 ( nQ ) 2 ~ are purely antiferro- 
magnetic, the orientation of n in the (three-dimensional) 
3-4-5 hyperplane determining the direction in real space 
of the antiferromagnetic (i.e., Neel) vector order param- 
eter. The novelty of Zhang's approach lies in its assem- 
bling of these two order parameters into a unified or- 
der parameter n, and the consequent possibility of ori- 
entations of n that do not lie wholely in one or other 
of the superconducting and antiferromagnetic subspaces, 
instead simultaneously containing components from both 
subspaces and, hence, characterizing regions that are at 
once partially superconducting and partially antiferro- 
magnetic. The purpose of the present Paper is to point 
out a simple but potentially interesting property of this 
model: in antiferromagnetic regions of the phase diagram 
this model supports three- (spatial)-dimensional antifer- 
romagnetic hedgehog configurations that find it energeti- 
cally favorable to have superconducting cores, as depicted 
schematically in Fig. |IJ 

It should be noted that the subject of the present Pa- 
per is, loosely speaking, conjugate to that of the recent 
one [0 in which it is shown that, within the S0(5) ap- 
proach, the cores of vortices in the superconducting order 
parameter should not be singular, the mechanism for the 
evasion of a singularity being escape from the two super- 
conducting dimensions into the three antiferromagnetic 
ones. 



II. ANTIFERROMAGNETIC HEDGEHOGS 
WITH SUPERCONDUCTING CORES 



What we mean by antiferromagnetic hedgehog config- 
urations with superconducting cores are energetically- 
stationary spatial configurations of the order parame- 
ter n(r) having the following properties: (i) Far from 
the (arbitrarily-located) center, the configuration n(r) 
closely resembles a (nonsuperconducting) antiferromag- 
netic hedgehog (i.e., a point defect in which the Neel 
vector points radially away from the center, or some 
global SO (3) rotation of this configuration) and the quan- 
tity sin 2 x [ = X)a=4( nQ ) 2 ]; which measures the de- 
gree of superconducting order (without regard to its 
phase), is small. Correspondingly, the complement cos 2 x 

[ = J2 a =i i n<1 ) }' w hich measures the degree of anti- 
ferromagnetic order without regard to its orientation, is 
close to unity, (ii) As the center of the configuration is 
approached, however, the order parameter escapes from 
dimensions 1, 2 and 3 into dimensions 4 and 5, so that 
superconducting order is acquired at the expense of an- 
tiferromagnetic order. Said equivalently, the angle x r0 ~ 
tates from to ir/2 as the center of the configuration is 
approached. (In principle, more exotic hedgehog excita- 
tions are possible, in which the antiferromagnetic order 
varies more rapidly. For the sake of simplicity we shall 
primarily focus on the simplest class.) By this mech- 
anism, the medium is able to remain nonsingular, and 
evade the (albeit finite) free-energy cost of the spatial 
gradient in the Neel vector (this gradient diverging as 
the center of the singular, purely antiferromagnetic, con- 
figuration is approached) at the expense of condensing 
locally into the "wrong" (i.e., superconducting) state, 
(iii) Whilst not being stable global-energetically — the 
homogeneous antiferromagnetic configuration of course 
having a lower free energy — antiferromagnetic hedgehog 
configurations with superconducting cores do turn out 
to be energetically favorable, compared with purely anti- 
ferromagnetic hedgehogs, as we shall see below, at least 
when amplitude variations of the order parameter are in- 
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hibited. Presumably, such configurations are also locally 
energetically stable ||. From the physical perspective, 
then, it would be quite intriguing if local regions of su- 
perconductivity were created by "stressing" the antifer- 
romagnetism in regions of the phase diagram in which the 
stable homogeneous state is not superconducting. More- 
over, the topological stability of these textures will tend 
to hold these "stresses" in place. 




FIG. 1. An eighth of an antiferromagnetic hedgehog hav- 
ing a superconducting core (determined numerically). The 
local orientation of the vectors indicates the local orientation 
of the antiferromagnetism. Their local magnitude indicates 
the local strength of the antiferromagnetism and, hence, the 
local weakness of the superconductivity. 

In the simplest version of Zhang's approach, the free 
energy F of a three-dimensional sample in which the or- 
der parameter n(r) varies with position r is given by 



F 



d A r 



u—l a—1 




(2.1) 



where p is the appropriate stiffness, v (= 1,2,3) runs 
through the cartesian spatial coordinates, and spatial 
anisotropies in the gradient term have been accommo- 
dated by coordinate rescalings. By tuning the chemi- 
cal potential /j, relative to its critical value p c (e.g., by 
doping), the parameter g [oc f/i 2 . — /i 2 )] is varied such 
that one moves from a region in which antiferromagnetic 
states are favored {g > 0) to a region in which supercon- 
ducting states are favored (g < 0) . This free energy is in- 
variant under separate rotations on the three-dimensional 
antiferromagnetic and two-dimensional superconducting 
subspaces; invariance under arbitrary five-dimensional 
rotations is absent whenever j / 0. Thus, from any 
configuration n(r) one can obtain a configuration having 
the same free energy via the transformation 



n -> (R A © R s ) • n, 



(2.2) 



where R A is a (3 x 3) orthogonal matrix operative in the 
antiferromagnetic (i.e., a — 1,2,3) sector (i.e., a magne- 
tization rotation) and R s is a (2 x 2) orthogonal matrix 
operative in the superconducting (i.e., a = 4,5) sector 
(i.e., a phase rotation), and the symbol © indicates that 
the five-dimensional operator is block-diagonal and com- 
posed of one three- and one two-dimensional block. 

To calculate the structure of an isolated antiferromag- 
netic hedgehog with a superconducting core, let us make 
the hypothesis that components of n(r) in this configu- 
ration can be expressed in the form 
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(2.3) 



Here, r, 9 and <f> are spherical polar spatial coordinates 
centered on the center of the configuration, and the func- 
tion x(r), which allows for interpolation between purely 
antiferromagnetic and purely superconducting values of 
the order parameter, is assumed to depend only on the 
radial distance from the center. This configuration is 
spherically symmetric, in the sense that for it we have 



n(i^ 



n(r), 



(2.4) 



where I s is the identity operation in the superconducting 
sector. By exchanging the radial variable r for the dimen- 
sionless version t (i.e., the radius, measured in units of 
the correlation length £ w = \f~p~Jg for the conversion of 
antiferromagnetic order into superconducting order) via 



X(t), 



(2.5a) 
(2.5b) 



and i nser ting the configuration (|]^) into the free en- 
ergy (2.1), we find that the free energy is given by 



- r (t 2 ■ t 2 

F = F j dt \ —X(t) 2 + cos 2 X(t) + — sin 2 X(t) 
Jo I 2 2 

(2.6) 

where F = 47rg(p/g) 3 / 2 , the overdot denotes a deriva- 
tive with respect to t, and y/p/gT is a large-distance 
cutoff, introduced to render finite the otherwise linearly- 
divergent free energy. Application of the calculus of vari- 
ations to the functional F then leads to the stationarity 
condition 



t 2 X + 2tX 



sin2AT = 0. 



(2.7) 
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FIG. 2. The degree of superconductivity sin 2 X as a func- 
tion of the scaled radial coordinate t (determined numeri- 
cally). 



The relevant solutions of Eq. ( |2.7| ) are: (i) X(t) = 
(i.e., the pure antiferromagnetic hedgehog, unescaped 
into the superconducting directions); and (ii) the solu- 
tion in which X(t) interpolates between ir/2 and as t 
varies from to oo (i.e., the antiferromagnetic hedgehog 
with a superconducting core). The precise form of the 
latter solution is readily found numerically, and is shown 
in Fig. ||. Its asymptotic behavior is {^-k — A) ~ t (for 
tC 1) and X ~ exp(— t) (for t — > oo). The configuration 
corresponding to solution (ii) is depicted in Fig. [lj. 

To determine which of the solutions, (i) or (ii), has 
the lower free energy, let us consider the quantity AF = 
_ F^), where and F^ respectively refer to 
the free energ y of solution (i) and of solution (ii) . Then, 
by using Eq. ( |2.7D , along with integration by parts, we 
find that 



AF = F 



X cos X + sin X 



| sin X, 



(2. 



where convergence at large t permits the replacement of 
the upper limit by oo. The numerical evaluation of this 
quantity gives AF ss — 0.272 F. This indicates that it 
is energetically favorable for the order parameter in the 
core of an antiferromagnetic hedgehog to escape into the 
superconducting directions. 



III. CONSEQUENCES OF AMPLITUDE-SECTOR 
FLUCTUATIONS 

As we have seen, in the setting of a model in which the 
amplitude of n(r) is constrained to be unity, there are 
hedgehog excitations that have superconducting cores. 
We now explore the issue of whether such excitations 
continue to exist in settings in which amplitude varia- 
tions of n(r) are inhibited (i.e., are not prohibited, al- 
though they are suppressed energetically). Under such 



circumstances, it is possible — and may prove energeti- 
cally favorable — for the core of the hedgehog to avoid 
antiferromagnetic gradient energy via the development 
of an amplitude-reduced purely-antiferromagnetic core, 
rather than by escaping into the superconducting direc- 
tions. To address this issue, we follow Arovas et al. (|] 
and consider a "soft-spin" generalization of the SO (5) 
model. Thus we consider a free energy of the form 



P / 



■ |(9 r n(r)) 2 + 2r 2 n(r) 2 cos 2 x(r) 

+ n(r) 2 (d rX (r)f + £" 2 n(r) 2 sin 2 X (r)} 

+ « / d\{- l -n{r) 2 + \n{rf}, (3.1) 

where 2a denotes the squared "mass" associated with the 
amplitude-sector fluctuations of n, n is the amplitude of 
n, and we have restricted the discussion to (spatially) 
spherically symmetric configurations. 

In the absence of amplitude fluctuations, the hedge- 
hog with superconducting core has n = 1 and \ varying 
from to 7r/2 as the center of the texture is approached. 
The amplitude-reduced purely antiferromagnetic hedge- 
hog excitation will be one for which n vanishes at the 
center of the texture and grows to unity at large dis- 
tances, and x = 0. The stationarity condition for n, 
which determines the structure of the amplitude-reduced 
hedgehog, can be solved numerically, allowing us to ob- 
tain the function n(r). Then, we may insert this back 
into Eq. (3.1) to obtain the free energy of the purely an- 
tiferromagnetic hedgehog with reduced-amplitude core. 
The free energies of the hedgehog with superconduct- 
ing core and the purely antiferromagnetic hedgehog with 
amplitude-reduced core each must be defined with a long 
distance cutoff to render them finite, but the difference 
between these two quantities is independent of this cut- 
off, and turns out to be given by 



Faf ~ F sc » 47T/9 (o.272y/p/g - lAtey 



(3.2) 



where the subscripts refer to the purely-antiferro- 
magnetic amplitude-reduced (AF) and superconducting- 
core (SC) hedgehogs. Thus, within the context of this 
model in which amplitude sector fluctuations are per- 
mitted, we find that the hedgehog with superconducting 
core will be energetically preferred when this quantity 
is positive, i.e., provided that g < 0.035a [or, equiv- 
alently, £ w > 5.35£ a , where £ a (= \J p/a) denotes the 
fluctuation correlation length for antiferromagnetic fluc- 
tuations] . Now, it is typical for £ a to be on the order of a 
lattice spacing for the cuprate materials, whereas £ w is ex- 
pected to grow as the superconducting phase boundary 
is approached from the antiferromagnetic state. Thus, 
one should anticipate that over a substantial portion of 
the antiferromagnetic part of the phase diagram, anti- 
ferromagnetic hedgehog excitations will have escaped- 
superconducting (rather than amplitude-reduced purely 
antiferromagnetic) cores. 
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IV. TOPOLOGICAL CLASSIFICATION OF 
HEDGEHOG EXCITATIONS 

We now turn to the issue of the topological classifi- 
cation of hedgehog excitations having superconducting 
cores, these excitations being nonsingular textures of 
the order-parameter field n(r). In pure antiferromag- 
nets, the existence of singular, hedgehog point-defect ex- 
citations is expressed, mathematically, by the statement 
n 2 (S l 2) = Z [|,||. What this means is that mappings 
(provided by order-parameter configurations) of spheres 
in real-space into the antiferromagnetic order-parameter 
space S*2 fall into homotopically inequivalent classes la- 
beled by the integers (and combine according to integer 
arithmetic). Within the SO(5) approach, however, the 
nonsingular hedgehog texture excitations having super- 
conducting cores are described by order parameter con- 
figurations n(r) that lie in the antiferromagnetic sub- 
space S*2 at large distances from the core, but escape 
into the full order-parameter space S4, as the core is 
approached. In order to complete the classification of 
these textures, then, we should ascertain whether or 
not there exist homotopically inequivalent textures that, 
at large distances, are homotopically equivalent and lie 
in the antiferromagnetic subspace Si- The appropriate 
mathematical machinery for this task involves relative 
homotopy groups and exact homotopy sequences 
To implement this machinery, we consider mappings of 
cubes (in real-space) such that the surface of the cube 
is mapped into £2 whereas the interior of the cube is 
mapped into S4. Such mappings are classified according 
to the relative homotopy group 113(54, 52). This group 
is readily computed by making use of the exact sequence 
of homomorphisms: 



n 3 (5 4 ) ^ n 3 (5 4 ,5 2 ) 3 n 2 (5 2 ) ^ n 2 (5 4 ) 



(4.1) 



Here, ^3, 73 and a 2 denote mappings of the elements of 
the previous group in the sequence to elements of the 
following group, that, in general, are not isomorphic 0. 
Now, as n^S^) and 112(54) are both the trivial group, 
the homomorphism 73 is, in fact, an isomorphism [js],^), 
from which it follows that 113(54, S2) — n 2 (S , 2) = Z and, 
thus, we find that there is no structure in 113(54, S2) be- 
yond what was already present in 112(52). The physical 
consequence of this result is that whilst hedgehog excita- 
tions fall into homotopically inequivalent classes, the pos- 
sible nonsingular superconducting cores of a given class 
of hedgehog are homotopically equivalent to one another. 



V. RELATED STRUCTURES IN OTHER 
CONDENSED STATES 

The notion of the conversion of singularities into tex- 
tures via the escaping of order-parameters into additional 
directions has been realized in several other condensed 
matter settings. For example, (uniaxially) nematic 



liquid-crystalline media have long been known to exhibit 
a structure closely related to antiferromagnetic hedgehog 
configurations with superconducting cores. When con- 
fined to a cylinder that imposes homeotropic (i.e., per- 
pendicular) boundary conditions on the nematic align- 
ment, the system can evade the threading of the cylinder 
by a singularity because the order parameter orientation 
can escape from the radial plane into the axial direc- 
tion H . This mechanism remains energetically favorable 
even for diamagnetic nematics in an axial magnetic field 
(for which escape costs condensation energy). 

Superfluid 3 He is another system that provides a rich 
array of topologically interesting textures |J. The ex- 
ample having the most relevance to the present Paper is 
that of hedgehog excitations in 3 He-B. The order param- 
eter for 3 He-B is a complex-valued 3x3 matrix of the 
form A^y = e'^-R^n, 0), where is a phase angle and 
i? M „ is a rotation matrix about the unit vector n by an 
angle 9. On long length-scales, 9 becomes fixed, due to 
a dipolar coupling, acquiring the value 9\ J1 known as the 
Leggett angle, so that the low-energy degrees of freedom 
are expressed by the possible values of <p and the direc- 
tions of n. Thus, the order-parameter space G is given 
by G = U(l) x 5 2 , so that n 2 (G) = n 2 (5 2 ) = Z, so 
that the system may form hedgehogs with the unit vec- 
tor n. On short length-scales, however, 9 can vary, so 
that the order-parameter space is effectively enlarged to 
U(l) x SO(3). We see that, as n 2 (SO(3)) = 0, over short 
distances there are no topologically stable point-defects, 
so that hedgehogs in 3 Hc-B have nonsingular cores. 

A similar effect occurs in nematic liquid crystals [ ^0[ , 
where on large length-scales the system is uniaxial, so 
that the relevant order-parameter space is RP2 (i-e., the 
real projective plane constructed by identifying oppo- 
site points on 52). On short length-scales, however, the 
order-parameter space is enlarged to 54, so that disclina- 
tions in the nematic order have nonsingular cores. 



VI. EXPERIMENTAL SIGNATURES OF 
HEDGEHOGS; CONCLUDING REMARKS 



We now briefly consider some issues associated with 
antiferromagnetic hedgehogs having superconducting 
cores that might be relevant to experiments. These ex- 
citations should be present after performing a quench, 
from high temperature or high magnetic field, into the 
antiferromagnetic state. The number of excitations per 
unit volume should be higher for more rapid quenches. 
As time proceeds after the quench, the number of hedge- 
hog excitations can decrease via their mutual annihila- 
tion, although this would require collisions of two or more 
hedgehogs. Presumably, this process can occur relatively 
slowly, at least at sufficiently low temperatures, so that 
one may anticipate regimes in which these excitations, 
once created, remain long enough for their consequences 
to be detected. 
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What experimental signatures might antiferromag- 
netic hedgehogs with superconducting cores yield? Let 
us suppose that a sufficiently high density of such excita- 
tions can be created, and that this density can be main- 
tained for a sufficiently long time. Then one may crudely 
regard the excitations as providing a set of randomly lo- 
cated, randomly phased, superconducting inclusions [ fll| . 
These inclusions would not be unlike Aslamazov-Larkin 
paraconducting fluctuations [Q , except that they would 
be "externally" maintained and, therefore, could be much 
longer lived. One might hope that these inclusions would 
be detectable in electrical conductivity experiments, their 
presence leading to an enhancement of the conductiv- 
ity. (One would need to account for scattering from the 
antiferromagnetic hedgehogs which, presumably, dimin- 
ishes the conductivity.) This enhancement should be sup- 
pressed by magnetic fields, and by the decay of the ex- 
citations. Similarly, one might also envisage observing 
Andreev reflection from the superconducting inclusions 
(although capacitive charging effects may suppress this 
effect ||). 

An externally applied magnetic field will be partially 
screened by these inclusions, leading to a negative con- 
tribution to the magnetic susceptibility. To estimate the 
size of this effect, we approximate the hedgehog cores to 
be uniformly superconducting and spherical in shape [ p"4"| . 
In the regime where the London penetration depth A is 
much longer than the core radius £, this leads to a dia- 
magnetic susceptibility contribution x = —S^/AOttX 2 , 
where S is the number of excitations per unit volume. 
One might also hope that the presence of antiferromag- 
netic hedgehogs with superconducting cores would be de- 
tectable via probes such as nuclear magnetic resonance, 
electromagnetic absorption and hedgehog/antihedgehog 
pair creation and, perhaps fancifully, scanning tunneling 
microscopy (e.g., with a magnetic tip). 

In addition, these excitations should leave their finger- 
print on the (staggered) magnetic structure factor S(k), 
this factor being determined by N(k) (i.e., the Fourier 
transform of the antiferromagnetic Neel vector at the 
probing wave vector k) via 



S(k) = V~ x N(k) • N(k), 



(6.1) 



where V is the volume of the system. Specifically, for 
length scales that are long compared with the core size 
£,r but short compared to the spacing between the hedge- 
hogs d^ 1 ^ 3 , we expect S{k) to have the conventional 
hedgehog form. However, for length scales that are short 
compared with the core size but long compared with the 
lattice spacing, we expect a reduction in S(k) (and hence 
scattering), owing to the diminution of the Fourier am- 
plitude of the antiferromagnetic moment at these scales. 
The computed hedgehog structure does indeed realize 
this scenario: 



However, it should be noted that, being sensitive only to 
antiferromagnetic order, this particular probe does not 
directly ascertain whether or not the cores of hedgehog 
excitations are superconducting (except via the depen- 
dence of the size of the cores on the location of the system 
in the phase diagram). 

If detected in experiments, antiferromagnetic hedge- 
hogs with superconducting cores would provide striking 
evidence in support of Zhang's SO (5) approach to the 
physics of high-temperature superconducting materials. 
Their presence would corroborate the notion that super- 
conducting excitations are essential low-energy excita- 
tions of the antiferromagnetic state. Moreover, it would 
prove rather intriguing to have at hand a physical sys- 
tem in which superconductivity is induced by the distor- 
tion of a thermodynamically preferred nonsuperconduct- 
ing state. 
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